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Abstract 

We compute the Green functions and correlator functions for N twist 
fields for branes at angles on T 2 and we show that there are N — 2 different 
configurations labeled by an integer M which is roughly associated with the 
number of obtuse angles of the configuration. In order to perform this com- 
putation we use a SL(2, M) invariant formulation and geometric constraints 
instead of Pochammer contours. In particular the M = lorM = N — 1 
amplitude can be expressed without using transcendental functions. We de- 
termine the amplitudes normalization from N — > N — 1 reduction without 
using the factorization into the untwisted sector. Both the amplitudes nor- 
malization and the OPE of two twist fields are unique (up to one constant) 
when the e <H> 1 — e symmetry is imposed. For consistency we find also an 
infinite number of relations among Lauricella hypergeometric functions. 
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1 Introduction and conclusions 



Since the beginning, D-branes have been very important in the formal de- 
velopment of string theory as well as in attempts to apply string theory to 
particle phenomenology and cosmology. However, the requirement of chiral- 
ity in any physically realistic model leads to a somewhat restricted number of 
possible D-brane set-ups. An important class are intersecting brane models 
where chiral fermions can arise at the intersection of two branes at angles. 
An important issue for these models is the computation of Yukawa couplings 
and flavour changing neutral currents. 

Besides the previous computations many other computations often involve 
correlators of twist fields and excited twist fields. It is therefore important 
and interesting in its own to be able to compute these correlators. As known 
in the literature pQ and explicitly shown in [2] for the case of magnetized 
branes these computations boil down to the knowledge of the Green function 
in presence of twist fields and of the correlators of the plain twist fields. 

In this technical paper we have analyzed the N twist fields amplitudes at 
tree level for open strings localized at D-branes intersections on T 2 using the 
classical path integral approach pp. The subject has been explored in many 
papers and in both the branes at angles setup and the magnetic branes setup 
see for example ([3J, [3|, [5], [6], [7], [8], [9]). 

We have shown that there are different sectors with different amplitudes. 
Sectors are labeled by an integer M cw (1 < M cw < N — 2) and that the 
number of sectors is equal to the number of obtuse angles formed by the 
brane configuration. This means that for example all the configurations in 
fig. ([TJ have different amplitudes. In particular the quantum amplitudes 
with M cw = 1 can be expressed using elementary functions only. This result 
generalizes the result previously obtained for both four point amplitudes (|6j, 
|7J) and for the N point amplitudes [5] where only the special case M = N — 2 
were considered. Since the N = 4 M = 2 amplitude has also been obtained 
by a different approach in [10]. it would be interesting to understand how 
this can come about in this different setup. 

We have also obtained the normalizations (up to one constant) of both two 
twist fields OPE and amplitudes. This result has been achieved using three 
ingredients: the consistency of the N twist fields correlator factorization into 
N — 1 twist fields one, the canonical normalization of the 2 twist correlator 
(a e (x)ai- e (y)) = 1/ (x—yY^^ and the assumption of the symmetry of under 
<r £ o (Tl_ e . 

Finally we have computed the Green functions in presence of N twist 
fields and we have shown that in order to do so there needs three different 
kinds of derivatives instead of the usual two which are needed in the closed 
string case. 

This paper is organized as follows. In section 2 we review the geometrical 
framework of branes at angles and we fix our conventions. In this section we 
discuss carefully how to make use of the doubling trick in presence of multiple 
cuts and the existence of local and global constraints. In section 3 we show the 
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Figure 1: The four different cases with N = Q. a) M ccw = 2 and M cw = 4 where 
M ccw is measured counterclockwise and M cw clockwise, b) M ccw = 3 and M cw = 3. 
c) M ccw = 4 and M cw = 2. d) M ccw = 5 and M cw = 1. 

existence of N — 2 different sectors and compute the corresponding classical 
solutions. We show also explicitly the results for the N = 3 and N = 4 
cases. Moreover using the known relation between closed string and open 
string amplitudes [TT] we express the classical action as a sum of products 
of holomorphic and antiholomorphic parts. Details on this computation are 
given in appendix [Aj In section 4 we compute the Green functions for the 
different sectors and give explicit expressions for N = 3 and N = 4 cases. In 
particular we discuss the existence of infinite relations among polynomial of 
Lauricella hypergeometric functions which must follow from the consistency 
of the procedure. Finally in section 5 we compute the quantum correlators 
of iV twists and their normalization factors. In particular we show that the 
M cw = 1 sector amplitudes can be expressed as a product of elementary 
functions. Moreover we discuss how N — 1 twist fields amplitudes can be 
obtained from N twist fields ones. A mathematically curious consequence is 
that certain determinant of order N — 2 involving Lauricella hypergeometric 
functions of order — 3 are expressible as product of powers. 

2 Review of branes at angles 

The Euclidean action for the string configuration is given by 

S = -— [ dr E I da (daX 1 ) 2 = -— *— [ d 2 u (dXdX + dXdX) (1) 
47ra' J J Ana' J H 
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where u G H, the upper half plane, d 2 u = drEda = du 2 f u and I = 1,2 so 
that X = -^{X 1 + iX 2 ), X = X*. The complex string coordinate is a map 
from the upper half plane to a closed polygon S in C, i.e. X : iJ — > £ C C. 
For example in fig. [2] we have pictured the interaction of N = 4 branes at 
angles Di with i = 1, . . . N. The interaction between brane Di and Di + \ is 
at fi G C where we use the rule that index i is defined modulo N. 

X(<t,t) 



a = 7r 



a = 




D 1 D 4 D 3 D 2 D 1 



n 

Figure 2: Map from the Minkowskian worldsheet to the target polygon E. 
2.1 The local description 

Locally at the interaction point fi the boundary conditions for the brane Di 
are given by 

Re(e- i7Ta *X' loc \ a=0 ) = Im(e- ma 'X loc \ a=0 ) - 9i = (2) 
while those for the brane -D^+i by 

Re(e- lwa ^Xi oc \ a=n ) = Im(e- i7Ta ^X loc \ a=n ) - g i+1 = (3) 

with 

ina i+1 _ i-irai 

fi = . 9 : 9 T (4) 

sin -K{a i+ i - (Xi) 

When we write the Minkowskian string expansion as X(a, r) = Xl(t + 
a) + Xr{t — a) the previous boundary conditions imply (and not become 
since they are not completely equivalent because of zero modes) 

X' L loc (0 = e^X' R loc (0, X' L /oc (£ + tt) = e i2 ™^X' R loc (£ - vr) (5) 
or in a more useful way in order to explicitly compute the mode expansion 

X' L loc {i + 2vr) = e^X' L loc (0, X' R loc (£ + 2tt) = e-^X' R loc (0 (6) 
where we have defined 

_ f ("i+l - «i) Oj+l > OLi , , 
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Figure 3: The connection between e and the geometrical angles as defining the 
branes. 

so that < €i < 1 and there is no ambiguity in the phase e l27T£i entering the 
boundary conditions. The quantity net is the angle between the two branes 
Di and Di + \ measured counterclockwise as shown in fig. [3] A consequence 
of this definition is that e becomes 1 — e when we flip the order of two branes. 
For example the angles in fig. [4] become those in fig. [5] when we reverse 
the order we count the branes, i.e. when we follow the boundary clockwise 
instead of counterclockwise the physics must obviously not change. 
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Figure 4: A polygon S with an ob- 
tuse angle and branes counted coun- 
terclockwise with N = 4 and M ccw = 
3. 



Figure 5: A polygon S with an ob- 
tuse angle and branes counted clock- 
wise with N = 4 and M cw = 1. 



We introduce as usual the Euclidean fields Xi oc (u, u), X[ oc (u,u) by a 
worldsheet Wick rotation in such a way they are defined on the upper half 
plane by u = e TE+ta G H. The previous choice of having brane Di at a = 
Q and brane -Di+i at a = ir ^ implies that in the local description where 
the interaction point is at x = Di is mapped into x > and Di + \ into 
x < 0. The boundary conditions Q can then immediately be written as 

dX loc (x + i0 + ) = e i2wa >dX loc (x - i0 + ) < x, 

8Xi oc (x + iO + ) = e i2nai+1 dXi oc {x-iO + ) x<0 (8) 
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and similarly relations for X which can be obtained by complex conjugation. 
When we add to the previous conditions the further constraints 

X loc (0,0) = f h X loc (0,0) = f: (9) 

we obtain a system of conditions which are equivalent to the original ones (|2j 

In order to express the boundary conditions (pi) in the Euclidean formu- 
lation it is better to introduce the local fields defined on the whole complex 
plane by the doubling trick as 

„ , . _ J dXi oc {u) z = u with 3z > or z G M + 
dX loc {z) - | ^ a ,Qx loc {u) z = u with < or z G M+ 

p.y / \ _ J dXi oc (u) z = u with 9z > or z G R + , . 

dX loc {z) - | e -i2^Q Xioc ^ z = u with Q Z < o or z G R+ ( ' 

In this way we can write eq.s Q as 

dX loc (^5) = e i2 ^dX loc (6), dX loc (e i27r 5) = e - i2 ^dX loc (5) (11) 

Notice that while the two Minkowskian boundary conditions ^ are one the 
complex conjugate of the other the previous Euclidean ones are independent 
and each is mapped into itself by complex conjugation therefore the Euclidean 
classical solutions for X and X are independent. 

The quantization of the string with given boundary conditions yields 

1 OO — qjT 

x l0C ( U ,u) = u +4v2^ v a{ ; )n u-("+i-«) - -<2lL„»+« 



2 *-~L n + 1 - €j n + ej 

?1=U 

1 00 Q! 

2 ^ n + l-ei n + a 

n=0 



,-,t 



*ioc(u,«) = /? +i-V2^V H „ 



2 ^n + l-ej n + ej 

n=0 

1 oo — q,T 

2 ^n + 1-^ n + e; vy 

?1=0 

with non trivial commutation relations (n, m > 0) 

[ Q (*)"' tt (i)«J = ( n + £ i)<W, [«(i)n, «(j) m ] = (n + 1 - £i)(5m,n (13) 

and vacuum defined in the usual way by 

U(i) n \Ti} = a(i)n\ T i) = n > (14) 

The vacuum is then generated from the twist operator a €i j t which depends 
both on the twist e, and on the position /j G C. The dependence on the twist 
can be read f.x. from the OPEs 

dX{u)a tiJi {x) ~ (« - x)^- 1 (aXa eiJ J(x) 

aX(n)a eii/i (x) ~ (« - x)- £ H5Xa ei!/i )(x) (15) 
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which can be deduced from the local computations 

dX loc (u)\Ti) ~ u^ 1 (-fiV2^4 )0 |25)), a^ oc (n)|T,) ~ iT ei (-^ V^c^jT;)) 

(16) 

On the other side the dependence on /j can be read from the OPE 

e ihX ^a £iJi (x) ~ Izr^^e-s^^^^e^ffe,,/,^) (17) 
which can be deduced from the local computation 

i i-|a'* a e -JJI a ( ei )a / * a . e ik-X loc (z,z) . | T .\ ^ \ z \-±a'k 2 e ik-fi e -iR 2 (ei)a'kh T \ 

(18) 

upon the identification [33] e iA>x ( 2 ^) o |,z|~2 a ' fc e ~4 R ( e iW k • e ik-x loc (z,z) . 

with R 2 {ei) = 2^(1) — V^i) ~~ ^(1 ~~ e «); V^- 2 ) = ^dz^ being the digamma 
function. Notice that there is no obvious way of computing the angles aj and 
from OPEs. 

2.2 Global description 

In the local description, where the interaction point is at x = 0, Di is mapped 
into x > and D- l+ i into x < this means that in the global description the 
world sheet interaction points are mapped on the boundary of the upper 
half plane so that Xj+i < Xj. The global equivalent of the local boundary 
conditions eq.s ^ become 

dX(x + i0 + ) = e i2wa *dX(x - i0+) x { < x < Xi -x 

dX(x + i0 + ) = e- i2 ^dX{x - i0 + ) Xi < x < x t -i (19) 

To the previous constraints one must also add 

Xi oc (xi,Xi) = fi, Xi oc (xi,Xi) = f* (20) 

in order to get a system of boundary conditions equivalent to the original ones 
((2j [3]). When we introduce the global fields defined on the whole complex 
plane by the doubling trick as 



dX(z) 



dX{u) z = u with 9z > or z G R - [xn, X\\ 
e i2lTOll dX{u) z = u with 9z < or z £ R - [x N , xi] 



„y, , _ J dX{u) z = u with > or z G K — [xjf, x\] , . 
W ~ \ e- i27Tai dX(u) z = u with Qz < or z G R - [x N , x t ] ^ ' 

the local boundary conditions ^ can be written in the global formulation as 

dX{xi + e i27T 5) = e i2lT€ 'dX{xi + 5) 

dX{ Xi + e i27T 6) = e- i2w£i dX( Xi + 5). (22) 
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For the proper definition of the global constraints which follow from eq.s (20), 
for example when dealing with the derivatives of the Green functions as in 
section [4] it is worth noticing the behavior of the previously introduced fields 
under complex conjugation when z is restricted to z G C — [—00, x\] 



\dX{z))* = e- i27rai 8X(z -> z) = dX(z) 



dX(u) z = u 
e -«*K<*idX(u) z = u 



[8X{z)]* = e~ i2 ™-dX{z -> z) = dX(z) = I e i2^^ {u) J = I (23) 

where dX{z — > z) means that the holomorphic dX(z) is evaluated at z. 
The previous expressions also show that it is not necessary to introduce the 
antiholomorphic fields dX(z) and dX{z) which it is possible to construct 
applying the doubling trick on dX(u) and dX(u) respectively. 



3 The path integral approach 

Following the by now classic method [T] we compute twists correlators by the 
path integral 

(a eiJl (x 1 )...a eNjN (x N ))= [ VXe- SE (24) 

where Ai({xi, e^, fi}) is the space of string configurations satisfying the bound- 



ary conditions (H9| and (20). Since the integral is quadratic we can then 



efficiently separate the classical fields from the quantum fluctuations as 

X(u,u)=X c i(u,u)+X q (u,u) (25) 

where X c i satisfies the previous boundary conditions while X q satisfies the 
same boundary conditions but with all fi = 0. After this splitting we obtain 

(a tlJl (si) . . . a eNjN (x N )) = M( Xi , ei ) e - 5 ^(^i,/i) (26 ) 



The explicit expressions for X c i and X c i given in eq.s (32, 33) show that 
they vanish when fi = hence also the classical action evaluated for fi = 
SE.d( x i, tii fi = 0) is zero. Actually because of translational invariance what 
said before works even when all fi are equal, i.e. when fi = f and therefore 
we can identify 

M(xi, Ci) = {(T ei ,h=f{xi) ■ ■ ■ Ve N ,f N =f{x N )) (27) 

Our strategy is therefore first to compute the classical contribution in the 
rest of this section and then compute the quantum contribution in section [5j 
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3.1 The classical solution 



We want now to write the general solution for dX and dX in a way that the 
SL(2, R) symmetry is manifest. To this purpose we introduce the anharmonic 
ratio 

(z - x N )(x 2 - Xl) 

[Z - X 1 )(X2 ~ X N ) 

and the corresponding ones Uj where z has been replaced with Xj. In par- 
ticular we get lon = 0, u)2 = 1 and uj\ = —00. The choice of uj\ = —00 is 
dictated by the request that powers are defined as (w — = \u> — ijJi\ e e % ^ e 
where <fi = arg(u — Wj) is counted from the real axis with range (— vr,7r) so 
that all cuts must be towards —00, as it is shown in fig. Q 




co 3 oo 2 



Figure 6: Cuts and prevertexes positions in the u> plane. 



We can now write the general solutions as 

N-M-2 



dX(z) 



dz 



n=0 



M-2 



r=0 



(29) 



where we have defined the basis 

N 



d UJ X^ n \oj z ) = Y[(u z - w i )~ (1-e, ' ) gj™, 0<n<N-M-2 



J'=2 
N 



d u X^ (oj z ) = Yl (u z - ojj)- e 3 u r z , < r < M - 2 
i=2 

and we have also defined the integer 



(30) 



N 



M = Y' 



(31) 



i=l 



When the target polygon S is followed counterclockwise this integer M is 
equal to the number of obtuse angles plus 2 since every acute angle internal 
the target polygon is ir — ire while every obtuse one is 2ir — ire as shown in 
fig. ([7]). In a similar way when the target polygon is followed clockwise M 
is the number of acute angles minus 2. Nevertheless it is important to notice 




7T6,- 




D; 



D, 



+1 



Di 



D 



i+1 



Figure 7: When the target polygon E (the shaded area) is followed counterclockwise 
keeping the interior on the left side an internal acute angle is equal to ix — ne while 
an obtuse one to 2tt — Tie. 

how polygons having M and M' = N — M both measured counterclockwise 
or clockwise are not the same polygons as it shown in fig. ([!]) in the case 
N = 6. To distinguish between these two cases it is necessary to compare 
expression (|7j) with the phases ctj as derived from the geometrical relations 



fi+i - fi = e 



I fi+l - fi\ as shown in fig. fe 




Figure 8: The connection between f i+ i — fi and the geometrical angle an+i defining 
the brane. 

Also when changing the fs while keeping fixed the es the shape may 
change as shown in fig. Q for iV = 4 and M ccw = 2 and in fig. (10) for 

= 4 and M ccw = 3. From now on we measure M clockwise in not otherwise 
stated. Since the number of obtuse angles must be less or equal than N — 3 
we deduce that 2 < M ccw < N — 1 or 1 < M cw < N — 2 and hence there are 
N — 2 different sectors^ 



1 The symmetry [X c i(u,u;{l - e}, {/*})]* = X c i(u, u; {e}, {/}) maps M ccw into M c 



N - 



Mccw because it is like the map X — > X* which reverses the order in which a circuit is followed. 



Hence it does not map M ccw into M' cc 



N - M r . 
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Figure 9: The four different cases with N = 4 and M ccw = 2 and M cw = 2 which 
can be obtained moving the brane whose intersection points are the empty circles. 




Figure 10: The four different cases with N = 4 and M ccw = 3 and M cw = 1 which 
can be obtained moving whose intersection points are the empty circles. 



ensures the proper transformation 
( 1-€ j) and the corresponding 



- ti^- 



In the previous expressions (29) 
under SL(2,M) and the product IX^C^ 
one for dX yields the proper monodromies around all the point, x\ included. 
The extrema of the summations, i.e. the maximum allowed values of n and 
r are chosen in order to have a finite action and in particular their values are 
determined by the analysis of the behavior of the solutions around z = x\ 
and not around z = oo as one would naively expect. This happens because 
the solutions (29) behave as O (^) at z = oo because of the factor 

(cj. 



The powers of the products 11^2 



-(1- 



') and Yl 



N 



dui z 
dz 



uij)- 



are chosen in order to get a finite X c i(u,u) at the singular points, explicitly 
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using the definitions (21) and the expressions (29) we can write^] 

duj ? 



N-M-2 u 
X d (ll,u)=f 1 + ^ a n{uj) j 

n=0 Jx ^ 



dz 



dz 



N 

n 

3=2 



UJ, 



M-2 



h 



r=0 

N-M-2 

n=0 
M-2 



x±;z£H~ 



a N 



C?2 



a;, 



a;,- 



i=2 



duo Y[( 



UJ — UJj 



j r yUlj 



r=0 



j=2 
N 

duo Y\ ( w ~~ w i 

i=2 



(32) 



where in the last step we have used the explicit definition of the power to 
connect the integral performed in lower half plane with that performed in the 
upper half plane. In a similar way we can write 



N-M-2 



n=0 



N 



M 

11 £ M- 



dw Yliu-Uj)-^ 1 -^^ 



M-2 



+e 



-i2na-i 



UJ; 



r=0 



J=2 



N 



UJ — UJj) 3 bJ 



i=2 



(33) 



where the coefficients are again a and b and not a* and b* as naively expected 
because we computed both X d and X d using the definitions of dX and dX 
(21) which mix both dX and dX . On the other side X c i = (X c i)* hence 
there are constraints on the coefficients a and b, i.e. a* = e~ t2irai a and 
similarly for b but these constraints are precisely the ones needed to solve 
the equations (34) when one takes into account the geometrical requirements 
that f i+1 -fi = e ™+i|/ i+1 _ / f | as shown in fig. (§. 

In order to determine the N — M — 1 functions a(uij) (j ^ 1,2, N) and 
the M — 1 6(cjj) we need simply to impose the N — 2 geometrical constraints 



X c i(x i+ i, x i+ i) - X d (xi, Xi) = f i+ x - fi i 



N - 1 



(34) 



There is actually one more equation one can obviously impose, the one with 
i = l but it turns out to be linearly dependent on the previous ones when the 
geometrical constraints on / and e are imposed. It is worth noticing that the 
previous constraints have an obvious geometrical meaning differently from 
the the use of Pochammer path used in the literature. The explicit solution 



Because of way we have chosen the cuts we have [(ui z — uij) a ]* = (^z — ^j) a when cjj is r 
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of the previous constraints is given by solving the linear system of N — 2 
equation^] 

N-M-2 M-2 

£ (-r 1 I^\l-e j )a n +J2lS\^)br = e i ^=^(fi-f i+ i) i = 2,...N 

n=0 r=0 

(35) 

where we have introduced the real valued integral^] 

i n r u * N 

4,n ) K)= / ni^-^-r^" (36) 

which are connected to type D Lauricella generalized hypergeometric function 
when i,i + 1 7^ 1 by 

N 

4,0 w = n n - 1^ - ^+ii i_ai+i 

1 i-1 / N -oca N , 



dtr°»i (i-r^TT 1- + 1) TT 1- + * 

N 1 

n|, . , , \—OLj I, . , , I 1— Q! 7_|_1 — OLj 

, . M1 -6(1 - CCj, 1 - «i+l) 



• f d (1 - a«+i; 2 - - «j+i;?a) (37) 

where the parameters £ a (a = 1, ... N — 3) are given by 



J 1 < a < i - 2 

i-l<a<iV-3 



t — ) ^a + l-^i + l — — /qn\ 

Sa — S cjj-^j+i i ^ „ at o V JC V 



In particular for iV = 3 is Euler Beta function B and for iV = 4 F^ is 
a plain hypergeometric 2 -Pi, explicitly the previous expressions become 

/g(l-e i )=5(e 2 ,£ 3 + n) 

- £j ) = u,r 1+n (1 - ^P^ 3 " 1 2^(63; 1 - n - e 4 ; 62 + e 3 ; ^) 

4> - = ^ +e4+n_1 p/ * v 2^1(1 - 6 2 ; n + e 4 ; 6 3 + 64 + n; a*) 

i3(64 + n,e 3 j 

(39) 

3 The net effect of using the real valued integrals 1^' is simply the appearance of the phase 
e i7r £j =2 £j in the right hand side. 

4 All these integrals can be expressed using iffl since ujn = 0, explicitly (otj) = j}^ {ctj — 
fij,N n ) hut we have introduced this redundancy for notational simplicity. 

Moreover we have used w_/v+i = u>i = —00 because indexes are defined mod N. 
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We are now ready to compute the classical action for our solution. Using 
the explicit expression for dX and dX we can write 



N—M—2 „ N 



Sal = ~ ^7 [ J2 ( e ~ inaia n) (e~ lwai a m ) f d 2 u f[ \lo - u^ 2 ^ oo n O m 

n,m=0 ^ C j=2 

M-2 . N 

+ (e i7Tai b r )(e i7Tai b s ) / d 2 w JJlw-^T 2 ^ w r w s ] (40) 

r,s=0 ^ C j=2 

where an overall factor ^ appears because we have extended the integration 
domain from the upper half plane to the whole complex plane. Notice that 
e~ l7Tai a, e tnai b G M which is however not enough to use their moduli \a\ and |6| 
in the previous expression. As explained in appendix [A] using the technique 
developed in [H] the previous integrals can be expressed as a product of 
holomorphic and antiholomorphic contour integrals as 



N 

/ d 2 LO Y[\u-Uj\~ 2€] U) n U m 



i=2 

N-l N ( I 

= EE-k ) iff (e)ig?(6) (41) 

i=2 l=i+l \ j=i+l 

3.2 The explicit N = 3, M cw = 1 (M ccu; = 2) case 

Let us start examining the M cw = 1 computation. In this case we see imme- 
diately that dX is identically zero and that the only unknown is ao which is 



not a function but simply a constant. The eq.s (48) reduce simply to 



a e im2 B(e 2 , e 3 ) = h - h (42) 

(3) 

because „(ay) = B(l — a 2 , 1+71 — 03) where B(e2, £3) is Euler beta function. 
The complete solution is then 

X^\u, u) = h + (h ~ h) ffi £3) , 2F 1 (e 3 , l-6 2 ; 63 + 1; oj u ) u« (43) 

B{e 2 , e 3 J 

Consider now the M ccw = 2 case where only bo is different from zero and 
therefore dX = 0. Proceeding as before we get 

b e- i2 ™ij m *B(l -€ 2 ,1- e 3 ) = h ~ h (44) 
from which follows 

X f 2 -\u, u) = f 3 +(f 2 -h) (1 ' 1 ~ £3) . [2^1(1-63, e 2 ; 2-e 3 ; w u ) ^ 

_«(1 — 62, 1 — £3j 

(45) 

which explicitly shows the equivalence [X^'^ N M ^ cw \u, u; {1 — e}, {/*})]* = 
X^ MaM \u, u;{e},{f}). 



5 i-KOL-\ 



a n is real as discussed before but it is by no means assured that it is positive. 
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3.3 The explicit N = 4, M cw = 1 (M ccw = 3) case 

In this case we again immediately realize that dX is identically zero and that 



the only unknowns are the two functions 00(^3) and a\{oJz). The eq.s (48) 
reduce simply to 



1 



EH" lj S(l " £ 3) anM = e^=^(fi - f i+ i) 1 = 2,3 (46) 



n=0 



with I^S given explicitly in eq.s (|39|) The classical solution then reads 



n=0 



(47) 



The X^' 3cc (u, u) solution can then be obtained as xj£ (u, u; {e}, {/}) = 

[X^ l \u,U- {1-6}, {/*})]*. 

3.4 The explicit N = 4, M = 2 case 

In this case M can be understood either as M ccw or as M cw which of the two 
can be only decided looking at the phases {ct^}. The unknowns are the two 



functions 00(^3) and ^0(^3) and eq.s (48) reduce simply to 



(-) l - 1 /f ) (l - 63) aoM + 4° ) (6 3 ) b r (u 3 ) = e^= 2 ^ - f l+1 ) i = 2,3 

(48) 

hence the classical solution reads 

poj u 

X^ 2 \u,u) = h +a (LU 3 ) / dw (co - l) 62 " 1 ^ - ^r^V 4 " 1 



du {u - l)" £2 (a; - Lo 3 )~ t3 LJ~ t4 



oo:u) 



(LH+ 



(49) 



3.5 Wrapping contributions 



The wrapping contributions have been studied in [Tl] for the N=3 case and 
in [5] for the case M = N — 2 and there is not any difference among the 
different M values therefore the results obtained there are valid. 



4 Green functions for N > 3 

Having determined the classical solution we now compute the Green func- 
tions in presence of twist fields both as an intermediate step toward the 
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computation of the quantum part of correlators and as a key ingredient to 
the computation of excited twist fields correlators. 

Following partially the literature we define the following quantities for the 
quantum fluctuations which are connected with the derivatives of the Green 
functions as 

r v \^ - {d^q{z)dX q [w)a ei j{x x ) . . . a €N j{x N )} 

9(N,M){ z i w i\ x il) — 



( cr euf(xi)...a eNtf (x N )) 
(dX q (z)dXg(w)a eiJ (x 1 ) . . . a €N j(x N )) 
(a eitf (x 1 )...a eN j(x N )) 



/ <y ,„. {9Xq(z)dXg(w)a euf (x 1 ) . . .a £N j(x N )) 

HN,M)K Z -, w i \ x if) — / 7 s 7 u K M 

K ' {^e 1 ,f{ x l)---^e N ,f( x N)) 

We do not need to consider functions involving antiholomorphic quantities 
because dX and dX are related to dX and dX as in eq.s (23). 

Quantum fluctuations are required to satisfy the boundary conditions 

Re(e- ina >d y X q \ y=0 ) = Im(e- inai X q \ y=0 ) = x l+1 < x < a*. (51) 

These conditions that can reformulated as a set of local constraints 

dX q ( Xi + e i2lT 5) = e i2 ^dX q (xi + 6), dX q ( Xi + e i27T 5) = e~ i27rei dX q { Xi + 5) 

(52) 

and as a set of global constraints 

X q {xi,Xi) = X q (x i+ i,x i+1 ), X g (xi,Xi) = X q (x i+1 ,x i+1 ), (53) 

In the spirit of what done in the previous section we use a SL(2,M) in- 
variant formulation and we write 

1 (w -w) 6 ' -1 N ~ M M 



j^l y w ■>' n=0 s=0 
J z diO w 

dz dw 



h {NM) {z,w-{ Xl }) = e- i2 ™^^- M^i) d u X^{ui z ) d^ s \u v 



r,s=0 
N-M-2 



dz dw 

n,m=0 



(54) 



where a ns (uj), b rs (tdj) and c nrn (uj) are unknown functions of the anharmonic 
ratios Wjyi, 2 ,jv- 

Let us rapidly review the ingredients of the previous construction. The 
factors ,K 2 , ^^jf and 9 £ z z d £™ are there to ensure the proper SL(2,R) 
transformations. The powers of the singular parts have been chosen in order 
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to reproduce the singularities of OPEs 



dX{z)dX{w) ~ -. + 0(1) (55) 

(z — w) z 

dX(z)a eJ {x) ~(z- x) e - 1 {dXa eJ ){x) 

dX(z)a €j (x) ~ (z - x)- e (dXa eJ )(x) (56) 



where (dXa e j)(x) and (dXa e j)(x) are excited twists. In particular eq.s (52) 



introduced to generate ( 52 ) constraints. 



are the same of eq.s (56) as it should be since twist operators have been 



The upper bounds of the summation ranges are fixed by request that 



singularities z — > x\ and w — > x\ are not worse than those in eq.s (56) while 
the lower bound is fixed from the z — > xn and uu — > xn limits. 

There is another consistency condition: when Xi — > Xj we must obtain the 
corresponding Green function with N — > N — 1. It is worth stressing that at 
first sight there is a further constraint. Usually the OPE between two twists 
is written as 



(xi - xj) e ^ M{ei,ej) a e . +ejifj (x j ) e, + ej < 1 
iWr l (xi-Xj)^ 1 -^ 1 -^ Af^ej) a ei+e .- 1)f . e t + e 3 > 1 

(57) 

with A4(ei,€j) = N(ei,£j) = 1. We will discuss that it is not possible to set 
both Ai and N to one in section I5T31 Now we would however comment on the 
fact that the previous expression is written without higher terms leading to 
the wrong impression that all the omitted terms are descendants. If it were 



true that the OPE (57) has no other primaries in the rhs this would imply 
that the derivatives of Green functions are analytic functions of the variables 
Xi too since the overall singularity in Xi — Xj due to the power factor would 
cancel between the numerator and the denominator. This is not true as the 
explicit computations show but in the M cw = 1, M ccw = N — 1 case and 
the reason is that the previous OPE involves actually an infinite number of 
primary fields with powers of OPE coefficients which do not differ by integers, 
explicitly for + tj < 1 

n, i . { Lv i )a, .(Lr,) - Xj)**** MU;.< j) a ei+ejJi {xj) 

+ Y,c k (xi - Xjr^^dXfa^j.^) + . . . 
k=i 

(58) 

where c& are certain numbers and . . . stands for other primaries and de- 
scendants. These primary fields have a simple interpretation as the states 
associated to the Hilbert space of twisted string since all of them have con- 
formal dimensions which differ by multiples of ±e. 

To continue and write in a more compact way the following expressions 
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we define 



P(oj z ,u w ) = J j 



N-M M 



S(oj z ,uj w )= y y a ns {uj)u^uj s w (59) 



n=0 s=0 



When we impose the constraint from the z —> w limit given in eq. (55) we 
get 



dS 



N-M M 


N 




n=0 s=0 


3=1 




N-M M 


s _ 1 - % 

3=2 


N 




n< 


n=0 s=0 


1=2 



or the equivalent equations with w — >■ z and 6^-1 — 9. These are (N + 
1) + N equations for a ns but only 2N are independent since both imply that 
djv-M, M = 0- Generically, i.e. for M cw / 1 these equations are not sufficient 
to fix the (N — M + 1)(M + 1) unknowns a ns and must be supplemented by 
the constraints which follow from eq.s ( |53| ). These further constraints allow 
also to fix the remaining (M — l) 2 + (N — M — l) 2 unknowns functions b rs 
and c nm . For example from the first equation in (53) we get 

X q (xi,Xi) - Xg(x i+1 ,x i+ i) = / dX q = I dx[dX q (x + i0 + ) + dX q {x - 

= f 1 dx[dX q (x + iO + ) + e i2 ™ 1 dX q (x-iO + )} = 

(61) 

which implies the constraint^] 

dx g( N ,M)( x + ^ + ,w) + e i2wai h iN)M) (x -i0 + ,w) = 

i 

dx l [NM) (z, x-iO+) + e i2nai g {NM) (z, x + i0 + ) = (62) 



Xi+l 



6 It is worth noticing that the segment [ajj+i, Xi] is followed for one addend above and for the 
other below the cut (it works also the other way round w.r.t. the main text). This ensures that 
both addends have the same phase modulus tt. Consistency among possible formulations of the 
constraints is due to [g(z,w)}* = g(z,w), [h(z,w)}* = e ii7Tai h(z, w) and [l(z,w)]* = e~ l ^ a H{z, w). 
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These can be explicitly written as 



M N-M . Ui duj N 

E w » E / n (u _l v II - "i)* -1 " n 

M-2 AI—2 . u . N 

+ j2 < e m^o / duj n> - u ^~ ej ^ = ° ( 63 ) 

s=0 r=0 Juj %+1 ;u>eH- j=2 

N-M M ^ ^ N 

e ^ e / u r, _ u) 2 - w i) _ej U ' 

n=0 s=0 Jui+iweH- {"z . =2 

N—M—2 N-M-2 N 



+ E ^ E C nm(Uj) P duf[(u> - Uj)*- 1 UJ m = 

(64) 



n=0 m=0 Ju t+ i;weH j=2 



As in the case for the classical solution only iV— 2 of the previous intervals give 
independent constraints, let us say i = 2, . . . N — 1. All these constraints are 
then sufficient to fix completely and uniquely all the coefficients. These con- 
straints are actually much more than needed since they equate a polynomial 
in co w or in u z to an analytic function. If we expand around u w = u z = oo 
and consider only the polynomial part we have enough equations to fix all 



the unknowns since to the previous 2N constraints in eq.s (60) we add M — 1 
equations in oo w times N — 2 intervals and N — M — 1 equations for uj z times 
N — 2 intervals. Actually the previous equations are already overdetermined 



in a since the 2N eq.s (60) and the [N — 2){M — 1) ones in uj w are sufficient 
for fix both a and b and similarly for the ones in u z therefore for consistency 
we suppose that this overdetermined system is consistent as well as all the 
remaining equations obtained from the polar part in uj w and uj z . As far as the 
consistency of the functions a determined in the two ways we have checked 
it in particular limits in the explicit cases treated afterward. Moreover all 
constraints derived from the polar part are polynomials in the integrals 



(36) since the functions a and b are solutions of a linear system whose co- 
efficients are precisely the I^ N 'S. Nevertheless it is easy to show that all 
constraints must be equivalent to a relation with polynomial coefficients in 



uij^i2N an d 1^ and at most linear in /( Ar+1 ) (see eq. (75)) with one of the 
parameters equal to 2. This can be seen as follows. It is possible to split 
9(n,m)( z > w 'i i x i}) i n a singular part and a regular one as 

9(N,M) ( Z > W > i X i}) = 9s(N,M) ( Z > W 'i i X i}) + 9r(N,M) W i i X i}) 

1 (u -uj-Yi- 1 N ~ M M 

9 S (N,M) (*, w; { Xi }) = _ H ' _ ^ Yl E a (o)-KK< 

^ ' jy^l ' w 3' n=0 s=0 

doj du N ~ M ~ 2M ~ 2 

9r(N,M)( Z ' W '{ x i}) = f, Z f, W E E UnsiUj) duX^ (u z ) d^X^ (uj w ) 



dz dw 

n=0 r=0 



(65) 



19 



This splitting is completely arbitrary and therefore it is not uniquely defined 
but it can be made unique imposing further conditions such for example the 
request of setting to zero all the CL n s= ] t — n but the twcpj with lowest n as 
for example in eq. (82) for the (N = 4,M = 2) case or the request that 
the singular part g s (N,M) S oes m t° 9s(N—l,M') when shown in 

appendix [B] and explicitly in eq. (150) for (N = 4, M = 2) case. Once fixed 
by a "gauge choice" the singular part the regular one is fixed by the global 
boundary conditions. 

Actually if we choose to split g into a regular part and a singular part 
(which is fixed not uniquely by the OPEs) as in eq.s (65) the equation (63) 
can be written as 



M 



s=0 



N-M 

a (0)ns(^j) 

n=0 
M-2 



2 



OJ - OJj ) 1 



N-M-2 



s=0 n=0 
M-2 M-2 



N 

n 

i=2 



Cj—1 , ,n 



UJ 



s=0 



r=0 



dojd u X {r \uj) = 



(66) 



which reveals that the singular part g s (N,M) contributes with a term linear 
in l( N+1 > (with one parameter equal to 2 because of the term (oj — uj w )~ 2 ) 
while the other terms have rational coefficient in uy^i 2,jv an d 1^' once we 
plug the solution for the coefficients back. In a similar way we can write the 



equation corresponding to (64) as 

N-M M 

E «*E 

a (0)ns\0Jj) 



dijj 



N 



71=0 s=0 

N-M-2 M-2 



(<*>z 



OJ 



Y[(u-u)j)- e i UJ n 



n=0 s=0 
N-M-2 N-M-2 
,m 



+ E 



OJ 



z 



brim (OJj ) 



J=2 

dojd^Xuj) 



m=0 



n=0 



eH- 



dood ul X^(u) = (67) 



Having determined the derivatives of the Green functions we can recon- 



7 Eq.s (60) divide naturally the set of unknowns {a ns } into subsets {a„ s } s -(-„ = fc and for each 
of these subsets there are two linear equations. 
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struct the actual Green functions as 



dv! 



dv' 9{n,m){u' ,v';{x})) 



+ I dv! / dv' l {N>M )(u',v';{x})) 

cv 



+ e 



dv! \ dv h( N ^ M )(u ,v'; {x})) 



Xj-.v'eH 



and 

r (N,M) 



Gf^ M) (u,u;v,v;{x}) 



+ / dv! / dv' g (NM) (u,v;{x})) (68) 

lxi]u£H~ JXj;v'£H 



u rv 

dv! I dv' 1( N ,m){u',v';{x})) 

Xi;u'£H Jxj\v'eH 
u fV 

i2n ai I du , 



+ e 



Xi;u£H JXj;v'£H 



dv' 9(N,M){u',v';{x})) 



+ r— i / dv' / dv' g {N , M )(v,u';{x})) 

J4na! I d -, 



Xi\udH 



Xj;v'£H- 



dv h {NM) (u ,v';{x})) 
(69) 



where the arbitrariness of the lower integration limit is due to the constraints 



(62). 



4.1 The explicit N = 3, M cw = 1 case 

In this case #(3,1) is completely fixed by the local constraints to be 

1 (u z - l) 62 " 1 



5(3,1) {z,w;{xi}) 



" (z - w) 2 (u w - l) e 2 

(1 - e 2 - e 3 )u;^ + (e 2 + £3)^0^ - (1 - e 3 )w 2 - e 3 w w (70) 



while /i(3,i) = since M = 1. We get therefore a constraint from eq. (63) or 
the equivalent form (66) which read 



.70 



a;,,, / <iw 



1 



+ / 



1 



(u - U} w ) 2 



{u>-lf*- 1 u e3 - L [(e 2 + e 3 )u-e 3 ] 

(00 - I) 62 " 1 ^[(l - e 2 - e 3 )^ 2 - (1 - e 3 )w] = 

(71) 
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which can either be read either as a constraint on the hypergeometric func- 
tions 



(l-e 2 -e 3 ) B(2 + e 3 ,e 2 ) 2 F 1 (2, 2 + e 3 ; 2 + e 2 + e 3 ; 



1 



(1 - e 3 ) B(l + e 3 , e 2 ) 2*1(2, 1 + e 3 ; 1 + e 2 + e 3 ; 



1 



(e 2 + e 3 ) B(l + e 3 , e 2 ) 2 Fi(2, 1 + e 3 ; 1 + e 2 + e 3 ; — ) - e 3 B(e 3 , e 2 ) 2 Fi(2, 1 + e 3 ; e 2 + e 3 ; — ) 



(72) 



or as infinite constraints on the coefficients of the uo w expansion which relate 
different Beta functions. We are now left to determine £( 3j i) from eq. (64), 
explicitly 



oo. 



+0J Z 



dto 



" {Uz-uf 



duo 



1 



+ / du 



o K - w) 2 
i 



(c;-l)- e2 w - e3 (l-e 2 -e 3 ) 



(w - 1)~ £2 W - £3 [(l " £2 - e 3 )w 2 - (1 - e 3 )w] 



(cj z - oo) 2 



:(u-l)- e2 ^ e3 (-e 3 ^) 

+c 00 f duo {oo- l) e2 ~ l w 63 " 1 = 
Jo 



(73) 



or 



c^(l-e 2 -e 3 ) /gj(e i; 2) 
+uo z [(1 -6 2 - 6 3 ) lg(e,;2) - (1 - e 3 ) I^%;2)] 
-e 3 i { 2 % r ,2)+c 00 4%,) = 



where we have introduced the function 

N 



a i (oo - uo w ) l3 uo T ' 

j=2 



(74) 



(75) 



which is a slight modification of our previous definition (36) and is still con- 

(n) 

nected to the Lauricella functions Fjj . In the uo z — > oo limit we can de- 
termine the unique unknown coefficient and hence the l( 3i i) normalization to 
be 



coo 



-(1 - e 2 - e 3 ) 



B(l - e 2 , 1 - es) 



(76) 



5(e 2 ,e 3 ) 

We get also infinite constraints from the subleading orders in uo z or plug- 
ging the previous value for coo back into eq. ( 73 ) an equation of the form 
Y, B 2 Fi = as eq. d72b. 



uo. b 



uo w = 
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4.2 The explicit N = 4, M = 1 case 



Again as the case before 5(4,1) is completely fixed by the local constraints 
only to be 



g(4,i)(z,w;{xi}) 



1 



£3 



J 2 



(2 - w) 2 (uj w - ( Ww - w 3 ) 

eiw^ + (1 - e\)u? z u w ~ [(1 - £3 - £4) + (1 - e 2 - e 4 )w 3 ]o;^ 

- [(e 3 + e 4 ) + (e 2 + e 4 )cj 3 ]u; z u; lu + (1 - e 4 )u;3u; 2 + e^u) w 

(77) 



and huu = since M = 1. As in the (3,1) case from eq. (63) or the 
equivalent form ( 66 ) we get the constraints 

- - £ 3 - 64) + (1 " 62 " 64)^3] 4,0 (! " £ i! 2 ) 

+ c^{(l - 0i) 4?(1 " 6i! 2) - [(e 3 + e 4 ) + (6 2 + e 4 )u;3] 1^(1 - e f , 2) + e 4W3 - e f ,2)} 



+ ei /g ) (l- ej ) + (l-6 4 ) W3 ^(l- ei )=0 

In particular notice that ~ i* 1 ^ is the Appell function. 

We can proceed to determine the £(4,1) function. This amounts to fixing 
the four functions cqq, coi, cio, en from eq. (64) which reads 

(-l) i+1 {^ 2 3 £l 4§(e i; 2) + u,, 2 (1 - ei)4?(e j; 2) - ^ [(1 - e 3 - e 4 ) + (1 - e 2 - ^^(l " ^ 2) 
- [(e 3 + e 4 ) + (e 2 + e 4 )w 3 ]4,i ) (e J -; 2) + u z (1 - e 4 ) w 3 ( e i5 2 ) + ^s^fife; 2 )} 



r( 4 ), 



(78) 



+ coo - ej) + coi Jf/(1 - e 3 -) + w 2 cio Ifo (1 - ej) + c u 2™(1 - e,) = 

(79) 

for i = 2, 3. When we consider the w 2 — )• 00 limit we get two sets of equations, 
the one from the coefficient of cj. 



r( 4 >/ 



r(4)/ 



(-i^/W^ + do^a-^-i-m /;;/( I -<,) -0 

and the other from the coefficient of co2 



r( 4 ), 



(-l) i+1 (l + ei ) /^(e,) + coo J$(l " <*) + coi - 6 



?(4) 



r(4). 







(80) 



(81) 



plus an infinite set of constraints from the coefficients of the polar expansion 
in ijj z or, equivalently plugging the previous value back into eq. ([79]) and 
equation of the form ( 2 iq) 2 F^ + = analogously to eq. (721). 



4.3 The explicit N = 4, M = 2 case 

This is the first case where here are more unknowns coefficients than equations 
from the local constraints and therefore we must use the global constraints to 
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fix completely 5(4,2) and determine both hu^) and 1(4,2) which are now both 
not vanishing. We can nevertheless fix the singular part g s u,2) by choosing 
020 = so we can get 



1 



3s(4,2) 



w 



l)e 2 



eiuo 2 z io w + (1 - e\)Lo z u? m - [(2 - e 3 - e 4 ) + (2 - e 2 - e 4 )w 3 ]u; 2 a; w 
+ [(1 - e 3 - e 4 ) + (1 - e 2 - e 4 )w 3 ]^ + (1 - e 4 )w 3 ^ + e 4 w 3 ^| (82) 



Using the global constraints for g and h as given in eq. (66) for i = 2,3 it 
is then possible to determine aoo (which corresponds to 020 after the split 
of 5( 4; 2) i 11 * a regular and singular part) and 600 ; i n particular taking the 
uj w — > 00 limit we get 



°oo 4fo(l - e ) - ft oo 45 ( e ) 



-(1-ci) /g(l-e 



[ aoo - 6) + 600 i$(e) = -(1 - ci) - r 

where the minus sign in the lhs of the first line is due a careful treatment of 
phases. In the limit u w — > 00 eq. (67) allows to fix coo and again aoo as 



[(1 - e 3 - e 4 ) + (1 - e 2 
[(1 - e 3 - e 4 ) + (1 - e 2 
(83) 



aoo I' 



$(e) - Coo J$(l - e) = -ei ig(e) 
a 00 4 4 ) ( e )+coo4 4 o ) (l- e ) = - e i4i ) (6) 



(84) 



The two previous ways of fixing aoo must be compatible and this can be easily 
verified at least in the w 3 — > 1~ limit. 



^)^ 3 ] 4 4 } (i 



e4 )u, 3 ] /$(!-€) 



5 The quantum twists correlators 

In this section we want to compute the N twists correlators in the N — 2 
different sectors determined by M. We can generically write the N twists 
correlators in the M sector as 

(11^ ,fM)) = fr ( X .- X A^ (85) 

i=1 lll<i<j<N\ x i x j) 

The powers A^ can be completely fixed as follows. From the proper behavior 
for Xj — > 00 we get the constraints = 2A(cr ei j.) = ej(l — e») where we 

have defined Aji = Ay for j > i. Now redefining A ^ A n3<j<j<Ar-i( c<; * ~~ 

w i) AiJ IlsKiKN-ii 1 - ^i)^ 21 ll3<j<Af-i u t iN an d remembering that W2 = 1, 
ujn = and (cjj — Wj) oc (x» — Xj) we can set all As to zero but Ah, A12, 
(3 < i < N— 1) and A2at which can now be fixed by the first set of constraints. 
Therefore we can choose a "gauge" where the previous correlator can be 
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written 

N 

(Jl^eiJiixi)) 



i=l 



n 



3<i<N 



(x 1 


- X2 )§[el(l- 




e»)+e2(l-e2)— ejv)] 






1 




(xi 


-xn) 1 ^ 1 




-£i)+ejv(l-ejv)-£2(l-e2)] 






1 





( X2 _ XAr )5[ e 2(l-£2)+ejv(l-ejv)-ei(l-ei)+Ei=3 1 ^C 1-6 *)] 
■ A (N,M)(Uj^l,2,N) (86) 

We can now proceed in the usual way. We first compute the expectation 
value of the energy-momentum tensor as 

{dX q {z)dX q (z)a eiJ (xi) . . . a eN j{x N )) 



«r(*)» 



then using the OPE 



lim g{z,w) - ^ 



(87) 



(88) 



we compute 



iV 



e,;(l - e, ; 



(;Z - Xj) 2 



(89) 



d Xi ln(Y[a ei; f(xi)) = Jim (z - sc* 
The function -A(jv,Af) i n the quantum case where fi = f can be determined 



from eq. (|85j) (j / 1, 2, iV) as 

AT 



i=i 



— xi dxj dijjj 



(90) 



The full expression for ^4(jv,Af); i- e - t ne case where all /j are generically 
different, can be obtained by multiplying the quantum value for the classical 



action as in eq. ( 26 ) . 



5.1 The M ccw = N - 1, M cw = 1 cases 



Using the expansion ( 59 ) for S and the constraints ( 60 ) we can easily deduce 
that 

2 



«no» 



1 / duj. 



2 V dz 



N 

E 

3=2 



N 



N 



(U Z -UJj) 2 



E^ II 



J=2 



Ul z — UJj 



3=2 



1 (Pi 

duj?„ 



UJ Z - UJj 



(91) 
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It then follows that (j / 1, 2, N) 

N 



§6jlog(n°*ei,/i=/0 C <)> 



i=l 



+ 



M-ei 



n 



1 duj d 2 S 



J 



oji dxj du> 2 



from which we can obtain A^ NfM ^ using eq. (90) to get 



dlnA 



(N,M) 



duj 



3 



n 



i & 2 s 



i 



2 4^.^'-^9a;2 



JV 



+ £ 



W M — OJ 3 — UJj 



UJj - LOl 



(92) 



(93) 



The main issue is then to compute 



in two cases, i.e. M r 



1 for which 



. This can be done immediately 
since the maximum ui w power 



is 1 and M ccw = N — 1 for which §-§ = since the maximum power is 1 



as it is obvious from eq. (59). In the former case we get 

AT 



i=l 



n 3 <i<iv-i^i-^)^ (1 " ei) 
n 3 <Kiv-i(i-^)~ e2e! 



(Xl — X2)^ £l( ' 1 ~ ei ' ) ~^-' li = 31 ei ( 1_ei )+ e2 ( 1 ~ <E2 )~ eiv ( 1 ~ eiv )] 



n 



3<i<Af-l w j 



( Xl _ X7v )|[ e l( 1 - e l)-E^ =3 1 e i( 1 - e i)+ e iv(l-<:iv)-e2(l-e2)] 



(x 2 — XN ^[^2(l-e2)+e N (l-e N )-e 1 (l-e 1 )+J2'iL 3 1 



(94) 



while in the latter we get the same result but with the substitution e — > 1 — e 
by expanding w z around The coefficients Cjv,i and Cn,n-i will be fixed 
in section 



5.4 and are given in eq.s (134). 



5.2 The JV > 4 and JV - 2 > M > 2 cases 



For all the other cases it is enough to use a slight modification of the technique 
used in [12] (see also [5]). 

The main idea of this approach is to define a new basis for the classical 
solutions (see eq.s (29)) and consequently for the non singular part of the 
derivative of the Green function g(z, w) (see eq. (65)) which are closed under 
certain operations needed to compute the correlators. 
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We start therefore by defining a new basis for the classical solutions 



N 



d w X^{ui) = H(lo - ojjp- 1 H (ui - Ul ) IeS 

3=2 /eSj 



N 



3=2 



Ul-Ulj) CJ 



\\ (w - uji) IeS 



(95) 



leSr 



where we have defined two ordered sets 

S = {N — M — 1 arbitrary different indexes chosen among 3, . . . N— 1} (96) 
and 

S = {M — 1 arbitrary different indexes chosen among 3, . . . N — 1} (97) 

and the subsets Si = S — {1} for any I e S and similarly for Sj. In order to 
be able to define the previous basis as a linear combination of the original one 
(30) we need that both n > and r > 0, i.e N — 2 > M > 2. In particular 
what follows works even if either Sj = or 5j = 0, i.e. M = N — 2 or M = 2 
for example when N = 4 and M = 2. 

We can now expand the regular part of g and h, I as 

dui z du v 

9r(N,M)( Z ' W >\ x i}) 



h(N,M)(z,w;{xi}) 



l(N,M)(z,W,{xi}) 



dz dw 

dui z dui u 
dz dw 

dui z dui u 
dz dw 



iaSIaS 



IeS Jas 



J^cjjfo) d u xW{w z ) d w x( J \u w ) (98) 



las J as 



Then we can find a solution of the first of the constraints in eq. (62) as 

N-2 



9(N,M)(z,W,{ x i}) = 9s(N,M)( Z > W >{ X i}) 
= 9s{N,M){z-,W,{Xi}) 



duj z 
dz 

8ui n 



i=l IeS J - 
JV-2 

»=i IeS 



dui 



. (ui - UJ W ) 2 

dui 
. {ui z - Ul) 



PS (ui,ui w ) 
^PS (ui z ,ui) 



(99) 



where Sq(ui z ,ui w ) is the same as in the second equation in (59) but for the 
singular part of g, i.e. with coefficients O(q) as i n ec L- (65). Moreover we have 
also defined the (N - 2) x (N - 2) matrix W 



Wi 



1 = I dui O ul X^(u + i0 + ) i = l,...N-2, IeS 



dui d^X^iui - i0 + ) i = l,...N-2, IeS (100) 



Again as in eq.s ( 62 1 it is important that the integration is once above and once below the 
cut as this ensures that both integrals have the same phase modulus it. 
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From this expression is immediate to compute the energy-momentum tensor 
expectation value which can be split into a singular part as in eq. (91) but 
with the substitution S — > Sq and a regular part as 



((Tr(w))) 



a N ~ 2 

dw 



2J U w ^) l i 



i=l /es 



duj 



UJ„ 



Uj) 



S (uj,uj w ) 
(101) 



where j g" Si means j G {2, . . . N} \ Si. If we consider J G S we can then 
evaluated 



lim (w - xj)((T r (w))) 



9uk 
dx., 



N-2 



i=i 



duj 



Xi (u - ujj) 3 tJ Uj^s^J - w i) 



S (uj,ujj) 
(102) 



Now following |12j we rewrite the integrand as 



Si -l 



-S (uj,ujj) 



Les 

(103) 



. The leading singularity is O ({uj — ujj)~ 2+ej ) because So(ujj,ujj) = as 
follows from the first equation in (60) when evaluated for uj w = ujj. Moreover 
when the left hand side is subtracted the leading singularity and multiplied 



for n 



N 

3=2 



UJ UJq 



v l-e, 



/ flies — ^j) we are w ith a rational function with 



poles at uji (I G S) whic h va nish at uj = oo as the right hand side. Because 
of the sum over i in eq. ( |l02[ ) the only Tj[ needed is 



<l-ej 



V- 

^ OJJ 



1 l i 

- un 2 ii ujj - 



n 



UJ] UJ] 



-dLS (ujj,ujj) 



(104) 



When we insert this value into eq. (102) and add the contribution from the 
singular part which has the same expression as eq. (92) with S — > Sq and 



9 There are two integrals which are actually divergent but their sum is however convergent. 
These integrals correspond to the intervals for which ujj is a bou ndary point. 

10 We restri ct to the case J £ S because otherwise eq. (102) would contain the sum over all 
I € S and eq. (103) would contain the sum over all possible d u .d w XW(ui) each with a non trivial 
coefficient. 
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e — > 1 — e since we have here d^ z Sq we get 

A? 



i=l 



(l-6j)(2- £j ) 
Xj - Xl 



2(1 - ej) 



1-e/ Af - N + 1 - ei 

E + 



xj - Xl 



XJ - Xl 



duj_ 
dxj 



N-2 



8=1 



(105) 



from which the dependence on Sq has disappeared but we are left with a 
dependence on d^jWf . Differently from what done in [|12] we cannot rely on 
fact that twists have both an holomorphic and antiholomorphic dependence 
in order to end the computation using 



N-2 



det W = 



i=i 



and 



wj — uji 



(106) 



(107) 



Instead we have to rely on the second of the (62) constraints (or better its 
complex conjugate which is has the same expression with the substitution 
zti0 + — > =pz0 + ). Analogously as before we require J G S and we get 



N 



d Xj log([[ <7ei,fi=f(xi)) 



i=l 



xj — X\ 



2ej 
dxj 



+ 



-M + ei 



XJ — Xl Xj — Xl 



E 

1^1,3 

N-2 



i=i 



leSj 



j 



1 



(108) 



then only HJ = J<ESnSwe can average the previous expressions (105) 
and (108) Into this average we can use the analogous expression of eq. (107) 



a UJ wi 



l-ej 



u> j — Ulj 



W/ - wj 



(109) 
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and 

N-2 



N-2 



J2 (w-'yjd^wf + £ {w-'Yjd^w- 



i=i 



i=i 



N-2 N-2 

dujdetw- E E^" 1 )^^ 1 - E E( w_1 )& 

ieSj i=l /esj i=i 



d^&eiW+Y, 77^77 + E 



U!j — LOl f-rf UJ j — UJi 



1110) 



to get finally 

a wj log i^M)^) = d Wi7 log (det W)~^ JJ (wj - w/) 5 II ( W J - ^ 



n ~ 



Ajj-lKl-e^Cl-eO+ejej] 
(111) 



11 



The previous equation is valid only for J G 5 n 5 but if, in either S or in S 
there is at least one further element than those contained in 5 n S or if 5 n S 
contains all the independent loj as in the N = 4 case, we can deduce that 

yl(jV,M)( w i) = const(det W)~? 



ord(I)<ord(J);I,J€S ord(T)<ord(J);I,JeS 

Yl (cjj - wjJ^-iia-^Xi-^+^Q] (112) 



2<j<l<N 

as a consequence of the independence of the result under a change of the 

elements of 5 and/or S. Under the change S S' = (S \ {Iq}) U {ii} the 

integrals W 7 ^^ i n ec i- s ( 100 ) transform as 



WJ - OJL ( ' 0J L - OJ Io W l 



VV (S')i - VV (S)i 
so that the transformation of the determinant 



;ii3) 



det W s , t s = det W s> g J\ 



Las,, 



W/ - LO L 



IH4) 



is what is needed to compensate the change Tl orc i(i)<ord(j)-i JeS' ~~ 00 j) 2 

W-ord^Kordijyj^&S^I ~ U j) 2 ■ 



11 



In the expression we have used ord(I) to indicate the order of / in the ordered set S. 
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The final expression for the N twists correlator in the M sector is then 



i=l 



n 3 <Kiv-i(i-^)^ [{1 - e2){1 - ei)+e2eil 



(xi 




ei)+ e 2(l-e2)-ejv(l-<:jv)] 




ll3<i<jv-i w i 


)(l-ejv)+£jejv] 




1 


-<5i)+ € iv( 1 - € iv)- e 2(l-<:2)] 


(aJ2 


— x^lMi-^+tA^l-ejv)- 





•(detW^-i [] (wj-wj)i n 

ord(I)<ord(J);I,JeS ard(I)<ord(J);I,JeS 

(115) 

Notice that the previous expression is true even if there is only one element 
in S in which case the product 13(^1 ~~ ^j) 5 ^ s simply 1. Similarly for the S 
5.4 we will fix the constant C(7v,M)( e )- 



u 



case. In subsection 



5.3 N — 1 amplitudes from N amplitudes 

We want to check the consistency of the results of the previous section. We 
do this by making x,+i coalesce with Xj and so deducing the N — 1 twists 
correlators from iV twists ones. 

We start noticing that from the (N, M) sector we can generically compute 
both (TV, M) = (N-1, M) and (N, M) = (N - 1, M - 1) sectors depending 
whether €j + €j + \ < 1 or £j + e J+ i > 1. Exceptions are the M = 1 case where 
only M = 1 is possible and M = N — 1 where only M = N — 1 = N — 2 is 
possible. 



5.3.1 (N, 1) into (N - 1, 1) case 



Starting from eq. (94) we can very easily take the limit xj+i — > xj. When 
we use 

ej(l - ej) + ej + i(l - ej+i) = ej(l - ej) + 2eje J+ i ej = ej + ej+i (116) 
and 

(x J+i - xj)(x N - x x ) x 2 - X\ . . 

wj-wj+i = t T7 r I- 1 - 1 ') 



(xj - Xi)(x J+ i - Xi) X 2 - XAT 



we find 

N 



N-X 



(Yl a ei,f( x i)) 



XJ+1->XJ 



i=l 



(xj-xj+i)-™ 1 -M(ej,e J+1 ) ( JJ ^,/(^)> 

(118) 



i=i 
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and the consistency relation for the normalizations 

C {N ,i)(e) = C(Ar-i,i)(e) M(ej,e J+1 ) (119) 

where e are the twists of the (N — 1,1) theory defined by ij = ej for j < J, 
ej = e,j + e j + i for j = J and ej = ej + \ for j > J. Actually all the previous 
equations work even when we consider the ujj — > oo (xj — >■ x\) limit. 

5.3.2 (N,N- 1) into (N - 1,N -2) case 

In a way completely analogous to that done in the previous subsection we get 

N N-l 
(H a ei , f ( Xi )) ~ XJ+1 ^ XJ (xj - x J+1 y^-^-^ Af(ej, ej+i) ( j[ a^fa)) 
i=i ~ i= i 

(120) 

and the consistency relation the consistency relation for the normalization 
coefficients 

C (JV)1) (e) = C (N _ ltl) (e) N{ej,€ J+1 ) (121) 

where e are the twists of the (N — 1,1) theory defined by ij = ej for j < J, 
ej = ej + ej + i — 1 for j = J and e~j = for j > J. Again all works in the 
ujj — > oo (xj — > x\) limit. 

5.3.3 (N, M) into (N - 1, M) with 2 < M < N - 2 case 



In this case we start from the general expression (115) and choose the sets S 
and S so that J G S, J S and J + 1 S then we now show that the new 
sets S and S are given by S = Sj = S \ { J} and S = S. 

The previous choices are dictated by the need of having a simple and 
clean way of computing the limit of det W s §. In particular while the interval 

vanishes dXg^ J ^ and dX^ become the new dX^^ and dX^ 

and dX( J ^ develops a not integrable singularity at ujj = ujj+i and gives the 
leading singularity of det W s §, explicitly we find 



det W5,fl~W$ x det W §>§ (122) 



with 



w $s)i=j-i ~ (wj - ^i) 1 ^™ 1 b (zji n - n - ^ 



(123) 



where !?(•, •) is Euler Beta function. Using these results into (115) with a not 
so short computation we find the expected result 

N N 

(II a tijM)) (xj - x J+1 )- ejtJ + 1 M(ej,e J+1 ) { ]J ^,/ I (^)) 

i=l i=l,ij^J 

(124) 
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and a relation among the amplitude normalizations and the OPE normaliza- 
tion in eq. (57) which up to a phase reads 



C(jV,M)(e) [B(ej,e J+ i)} 5 = C {N _ 1M) {e) M{ej, e J+1 ) (125) 

where e are the twists of the (N — 1, M) theory, i.e. £j = ej for j < J, 
ej = ej + ej + i for j = J and = e J+ i for j > J. 



It is worth noticing that the previous result (124) shows that eq. (115) is 
valid even when S has only one element. If we perform the reduction from 
this case, i.e. with (N + 1,N — 1) and we compare with the expression for 
the (iV, N — 1) amplitudes we deduce that 

det ^ ( 0,5) n (cjj-uj)- 1 * n (uj-un)^- 1 (126) 

ord(I)<ord(J);I,JeS 2<j<l<N 

where ^ is simply the matrix || Wf ||. In other words certain determi- 
nants of order N — 2 (card(S) = M—l = N — 2) of Lauricella hypergeometric 
functions of order N — 3 (since all Wl can be expressed using are a 

product of powers. This could point to that also the general det W s § may 
be expressed as an elementary function. 

For the special case where both S and S have just one element, i.e. for 
iV = 4, M = 2 a direct and little different computation is needed but the 
result is the same. 

For checking the consistency of the approach and of the normalization 
coefficients we determine in the next section it is worth considering the ojj — > 
oo limit. The result for the normalization coefficients in this case is based on 
the relation 

~B(ej, 1 - £1 - tj) „y(»-»'-')-« ietWsS (127) 

and reads 

C ( iV,M)(e) [B{ej,l-e l -Ej)}- 1 *=C {N _ l>M) {e) M{ej,ei) (128) 

with the new twists given by ei = ei + ej, ij = ej for 1 < j < J and ej = e^+i 
for j > J. 

5.3.4 (AT, M) into (N - 1, M - 1) case 

In this case we can choose the sets S and 5 so that J S 5, J S and J+l 5 
then it is possible to show as in the previous case that the new sets S and S 
are given by S = S and S = Sj = S \ {J}. 

In particular it is possible to find analogously as before that the determi- 
nant behaves in the xj — > xj +1 limit as 

det Wsj-wMj^^ det W SfS (129) 



33 



and the amplitude reduction gives 



N N 

<n^ 1 / 4 (^)>~(aj-x J+1 )-c i -^ i -^)< n ( i3 °) 



It follows a relation among the amplitude normalizations and the OPE nor- 



malization in eq. (57) which up to a phase reads 

C ( iV,M)(e) [B(l-ej,l-ej+ 1 )]-?=C( N - ltM _ 1) (e)Ar(ej,e J+1 ) (131) 

where e are the twists of the (N — 1,M — 1) theory, i.e. ej = ej for j < J, 
ej = ej + ej+i - 1 for j = J and ij = e j+ i for j > J. 

As in the previous subsection starting from the (N + 1,2) amplitude and 
reducing it to (N, 1) we deduce that 



detW, 



(5,0) 



n 



-1 



ord(I)<ord(J);I,J&S 



2<j<l<N 



(132) 



where W^ s ^ is simply the matrix || W/ ||. 

5.4 Amplitudes and OPEs normalization 

We normalize the 2-point amplitude as 

1 



(cr e (a;)(Ti_ e (y)) 



(x- j/)^ 1 - 6 ) 



(133) 



This normalization is not unique since any redefinition as a e —> lZ(e)a e with 
TZ(e) 1Z(1 — e) = 1 would work. In particular this kind of redefinition can 
only be seen in amplitudes with at least three twist fields since it leaves 
unchanged amplitudes involving two twist fields and an arbitrary number of 
untwisted fields therefore it cannot be fixed factorizing a 4 twists into an 
untwisted channel. If we require the normalizations to be invariant under the 
symmetry eol-e then all the normalizations are completely fixed (up one 
constant k and phases) to be 



C(AT,i) = k 



N-2 



N 



1/4 



c, 



(N,M) 



a 



(N,N-1) 



N-2 



JV-2 



nf =2 r( gj )r(i- ej ) 
r( ei )r(i- ei ) 

1/4 

n m 



1/4 



2 < M < N-2 



(134) 
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along with the OPE normalizations 



M{a,P) = k 
= k 

M(a,P) = k 
= k 



r(i - a) r(i - p) r(a + /3) 

T(a) r(/3) T(l-a-/3) 

r(i-a)r(i-/3)r(i- 7 )i 1 / 4 



1/4 



r(a) r(/3) r( 7 ) 

r(a) r(2-a-/?)l 1/4 



a + /3 + 7 = 1 



r(i - a) T(l - /3) T(a + /3 - 1) 



r(a) r(/3) rffl 
r(i - a) r(i - p) r(i - <$) 



i/i 



a + fi + 5 = 2 (135) 



which also respect the symmetry e h 1 - e as AA(a, /3) = ,M(1 — a, 1 — /?). 
It is at first sight surprising that there is not symmetry among the twist 
operators in the M ^ 1, N — 1 case but this is due to two reasons. The first 
is our choice of using a SL(2, M) invariant formalism which singles out some 
points and the second is that not all twist operators are on the same footing 
since some couples of twists sum to a quantity less than one while others to 
one bigger than one. These normalization are the "square root" of the ones 
found in [15] for the N = 4 closed string case and matches those obtained for 
N = 3 in the magnetic brane case in [9] and for N = 4 case in |10j . 

Let us see how we can get the previous results by exploiting the con- 



sequences of equations of the previous subsections such as eq.s (125) and 



(131) 



First we notice that we can always normalize the 2-points correlator 
as chosen because the generic normalization factor C(2 j i)(e, 1 — e) is symmet- 
ric in the exchange e o 1 — e hence we can redefine the twist operators as 

£ e /y / 'C , (2,i)(e,l-e). 

From the reduction (JV = 3, M = 1) to (N = 2, M = 1) with the help of 
eq. (119) we find that M(a,(3) = C( 3il )(a, (3, 7) with a + j3 + 7 = 1 has the 



cr, 



following basic symmetries 



M(a,P) =M(P,a)=M(a,l 



a 



P) 



and all the others which follow from them. 

In a similar way from the (N = 3, M = 2) to (N = 2, M 



(136) 



1) reduction 



and from eq. (121) we find 

N(a,p) =M(P,a) 



N(a,2 



a 



(137) 



Now we can consider the (N = 4, M = 2) to (N = 3, M = 1) reduction 
in two different ways. Either with (ei, €2, £3, 64) — > (ei, 62, £3 + €4 — 1) which 
implies 

C(4,2)(ei,e2,e3,e4)[-B(l - e 3 , 1 - e 4 )-B(e 2 ,e 3 + e 4 - 1)]2 = 7V(e 3 , e 4 )C( 3] i)(ei, e 2 , e 3 + e 4 - 1) 

(138) 
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or with (ei, £2, £3, £4) — > (ex, £2 + £3, £4 — 1) which implies 

C(4,2)(ei)e2,e 3 ,e 4 )[B(l - £3, 1 - £ 2 )-B(£4,£3 + £2 - I)] 5 = Af(e 3 , £2)C(3,i)(ei, e 2 + e 3 - 1, £4) 

(139) 

Now taking the ratio of the two previous equations and using the symmetries 
of M and J\f we are led to the minimal ansatz 

M(a, p) = k[T(a)T{l3)T{\ -a- P)] a [T(l - a)T(l - j3)Y{a + (3)} b 

Af(a,/3) = k[T(a)T(P)T(2 -a- (3)} C [T(1 - a)T{l - (3)T(a + /3 - l)] d 

(140) 

which gives an overconstrained system when plugged back into the ratio con- 
straint whose solution is a = —b and c = — d = ^ + a. This solution imme- 
diately yields both C( 3)2 ) and C(4 )2 ) ■ Imposing the symmetry e o 1 — e then 
selects a = —\. It is then easy to generalize to the full expressions. These 
can be checked in different limits also when we consider cjj —> 00 using to eq. 
(128b. 



A Details on rewriting the classical action. 

We want to give some details on the use of KLT technique for reducing the 
integral 

/+00 r+00 N 

dx I dyY[(x + iy-uj j ) a i +n i(x-iy-uj j fi +fi i 
-00 J— 00 

(141) 

with rij ,rij G Z to the sum of products of an holomorphic and antiholomor- 
phic integral. First we interpret the previous integral in y as a path integral 
in the complex plane Y = t + iy. In the variable Y the integrand has cuts in 
±(ujj — x), the main issue is then to properly define the phase of 

(x + Y - u) a {x -Y -uf = \x + Y - oj\ a \x - Y - ufe***®. (142) 



The proper choice is shown in fig. ( 11 ) and is constrained by the request that 



when Y = iy and a = a then cf> + <p = 0. 

We can then rotate clockwise the path in Y plane, change variables as 
£ = x + t, rj = x — t and then we can rewrite the integral as 



• p+00 r+00 N 
j( N \a + n,a + n) = - - df / drj JJ |f - Uj\ a i \q - u $ \ 

^ J — OO J — OO ■_r > 



J 

X e - i7Ta 3 6(ri-uj) e -iiraij 0(£-Mj) 6(uij-v) (143) 
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Y = t + iy 



Figure 11: Proper definition of angles <fi and <fi and therefore of phases when x — u > 
0. When x — lo < we substitute (x — to) — > —(x — oS). 

If a = a then we can proceed as in KLT. We fix £ and we exam the 77 integral. 
Each factor of the integrand can then be rewritten as 

+ e{-Uj + e +i ™> [ Uj - (r? - i0 + )p' (144) 



when we choose the phase in the complex 77 plane as in fig. (12), obviously 
other choices would do the job as well. In words this means that when £ < ojj 
we run above the cut from —00 to coj in the complex r\ plane while we run 
below the cut when uja < £. Hence the original integral can be written as 





Imr\ 











Figure 12: Definition of phase in rj plane in the range (— 27r, 0). 



J (7V) (o + n,a + 



n) 



2 ru, N 
=JV— 1 u ' i + 1 3=2 



. iV 

e' E '- ai / dr, Y[ {uij -vT j {ri - WjP (145) 

j=2 



where the path d is given in fig. (13). We can then deform the C\ path in 
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Figure 13: The path Ci in the complex r\ plane for < £ < u>i. 



order to integrate between — oo and uJi+i above and below the cuts. The final 
result is then 



J iN) (a + n,a + 



n 



N-l N I I 

- E E sin r E 1 

i=2 l=i+l \ j=i+l 

j = 2 



(146) 



B Fixing the singular part of g(z, w) in a 
consistent way with N — > N — 1 reduction 

Let us suppose that all coefficients c ns (ujj) depend on ojj (3 < j < N — 1) 
in an analytic way. We want to show that it is then possible to fix then in 
a recursive way starting from those of the N = 3, M = 1 case. This can be 
done considering two limits Xj — ^ xjy, i.e. Uj — > and Xj — > x\, i.e. ujj — > oo. 
Combining the two cases when ei + €j < 1 and ej + < 1 we get 



JAT-1,M)/. - 



"n— l.s 



_(iV-l,M) 



(o),e)a;j 



(N,M) 
c 0,s 

SN,M) 
-N-M,s 



(w,e) 



,(JV-l,Af) 
(AT-1,M) 



(w,e) = %_ M '_i s (^,e) 



when l<n<iV — M — 1, < s < M and where we have defined 

h = e% i = l,... j-1 

h = e-t+i i = j, . . . N — 2 

k eAr_i = e w + £j - 9{e N + ej > 1) 



(147) 



(148) 



and 



ei = ei + ej - 0(ei + ej > 1) 

e? = e t i = 2, ...j-1 

h = e i+ i % = j, . . . N - 1 



(149) 
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and similar relations between Co with cu and Cj with u. For example applying 
the previous formula to the N = 4, M = 2 case we get 

+ (1 - ei - e 2 )u^ - [e 3 + e 4 + (ei + e 3 )uj 3 ]uj z uj w + (1 - e 2 

+ (1 - €4,)UJ 3 UJ Z + e±U) 3 U) w \ 
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